function R = ScatteringRatePower(p,x)
    s0 = x.*p(2)./p(4); % I/I0, p(2) is a scaling factor between measured cooling power and cooling light intensity.
    % p(4) = I0
    Rn = s0.*pi*p(5); % p(5) = \Gamma (linewidth)
    Rd = 1 + s0 + 4*(p(3)./p(5)).^2; % p(3) = \Delta (detuning)
    R = p(1)*Rn./Rd;
end